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Abstract 

We have investigated the contribution of connected triple and quadruple excita- 
tions to the binding in Be 3 by comparing various coupled-cluster (CC) and truncated 
configuration-interaction (Cl) treatments with multireference Cl (MRCI) and full 
Cl (FCI) calculations. The CC method with single and double excitations (CCSD) 
produces results that differ substantially from more elaborate treatments, but most 
extensions to CCSD that account approximately for connected triple excitations 
perform very well. In contrast, good agreement with FCI for Be 2 can be achieved 
only with the highest level CC and MRCI methods. 

I. Introduction 

For the diatomic molecule Be 2 , the CCSD method gives a poor potential en- 
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ergy curve, but after inclusion of connected triple excitations (CCSDT) a more 
reasonable curve is obtained [1]. Tbe curve is further improved after consideration 
of connected quadruple excitations [2]. While the bonding in the trimer and higher 
oligomers of Be is quite different from that in Be 2 , Lee et qL [3] found in a recent 
study of structure and bonding in Be and Mg clusters that the CCSD method also 
did not satisfactorily reproduce the binding energies determined with large-scale 
multireference Cl (MRCI) wave functions, even though the CCSD geometries were 
fairly good. In further investigations [4,5], however, it was observed that agreement 
with MRCI was greatly improved when the CCSD method was corrected approxi- 
mately for the effect of connected triple excitations [6], using the CCSD+T(CCSD) 
method [7] or an improvement denoted CCSD(T) [8]. While both methods include, 
noniteratively, all fourth-order contributions of connected triple excitations, the 
latter adds one fifth-order contribution which is normally of the opposite (i.e. posi- 
tive) sign. As such CC treatments are considerably cheaper than MRCI, they could 
provide a cost-effective means of studying larger Be and Mg clusters. However, 
these methods fail for Be 2 [2]. Hence, it would be desirable to estabhsh that the 
apparently good performance of CCSD+T(CCSD) and CCSD(T) for the trimers 
and tetramers does not arise from some fortuitous cancellation of errors in the cor- 
relation treatment. Such errors could be due to the failure to treat the connected 
triple excitation contributions iteratively, to the neglect of many terms in the full 
CCSDT equations (including several fifth-order terms), or to the neglect of higher 
than connected triple excitations, particularly the connected quadruple excitations 
that first contribute to the correlation energy in fifth-order. Such errors are a con- 
ceivable problem in Be clusters because of the known multireference character of 
these systems caused by the quasi-degeneracy of the Be 2s and 2 p orbitals. 

In the present work we examine a number of different correlation treatments 
for the Be 3 cluster. We perform FCI calculations to provide an absolute comparison 
standard for all of the correlation treatments. We also perform CC calculations with 
full iterative inclusion of triple excitations (CCSDT) [9,10], several approximations 
to the full CCSDT treatment as proposed by Bartlett and co-workers [6,7,9,11] (as 
well as the CCSD(T) method [8]), and two noniterative CCSD corrected schemes, 
which are correct through fifth-order and account approximately for the effects of 
connected quadruple excitations [2]. Effects of connected quadruple excitations are 
also estimated with the CCSDT+Q and CCSDT+Q* methods [12]. We compare 
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these results with those of single reference and multireference Cl calculations. Brief 
details of the computational methods are given in the next section, some details 
being deferred to section III, which contains our results and discussion. Conclusions 
are given in section IV. 

II. Computational methods. 

Two Be atomic natural orbital (ANO) basis sets [13] are used in this work. 
The larger is the [4s 2 p Id] basis given in Ref. 3, although in the present work the 
3s cont ami nant function is included in the basis. The smaller basis is obtained by 
deleting the d ANO to give a [4s 2 p] set. 

The [4s 2p] basis has been used in conjunction with a variety of correlation 
treatments: FCI, MRCI, CISD, CISDT, CISDTQ, MBPT(4), CCSD, CCSDT, and 
various approximate CC schemes for handling connected triple and quadruple exci- 
tations. The particular CC schemes are described in the next section. FCI calcula- 
tions in the [4s 2p Id] basis are beyond our computational capabilities, but all other 
methods were applied in this basis. In all of our single-reference calculations (in- 
cluding FCI) the MOs were obtained from an SCF calculation; the highest-lying a\ 
and e' MOs (these would correspond to Is complement functions in a segmented 
contracted basis) were deleted. The MRCI calculations were based on a CASSCF 
wave function with all Be 2s and 2p orbitals in the active space: all CASSCF config- 
urations were used as reference configurations. In the CASSCF/MRCI calculations 
the Be Is orbitals were taken from an SCF calculation and frozen, and the high- 
est a'j and e 1 MOs were deleted after the SCF calculation. The Be Is electrons were 
not correlated in any of the calculations. 

Calculations were performed around the equilibrium bond length assuming a 
D^h. geometry. A rather large grid spacing (0.3 cto ) was used because the Be 3 poten- 
tial is very flat. Optimum bond lengths were determined by fitting the computed 
energy points to a quadratic form in 1/t. 

All MBPT and CC methods are rigorously size-extensive and for this example 
separate correctly to closed-shell atoms. Hence, at the separated atom limit, three 
times the energy of one Be atom is obtained. Moreover, since the core electrons are 
not correlated, all of the CC methods used in this work give the FCI energy at the 
separated atom limit. In the [4s 2p] basis set, the SCF, MBPT(4), and CC atomic 
energies are -14.572920, -14.612790, and -14.617580 E h . In the [4s 2p Id] basis the 
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atomic energies for these methods are —14.572942, —14.614290, and —14.618643 Eh 
respectively. MRCI, however, is not rigorously size-extensive. Consequently, for 
this method binding energies were computed using a supermolecule approach for 
the dissociation limit. 

The SCF/CI and CASSCF/MRCI calculations were performed using the 
MOLECULE-SWEDEN [14] program system. Some CCSD and CCSD(T) calcu- 
lations were performed with VCCSD, a vectorized CC program [15] interfaced to 
MOLECULE-SWEDEN. Other CC calculations were performed using ACES [16]. 
The FCI calculations were performed using a modified version of the Knowles- 
Handy [17] program. Calculations were performed on the NASA Ames Central 
Computer Facility CRAY Y-MP/832, the NAS Facility CRAY Y-MP/8128, and 
the Quantum Theory Project FPS-164/MAX. 

III. Results and Discussion. 

The results of calculations in the [4s 2p] basis are summarized in Table I. 
Energies are given explicity and differences relative to FCI are given in parentheses. 
One unfortunate consequence of using a basis set this small is that Be 3 is not bound 
at any level of calculation. As Be 3 is actually bound by at least 24 kcal/mol [3], 
this indicates a very significant coupling of the one-particle and n-particle spaces, at 
least for the binding energy. However, there is a local minimum in the energy at a 
bond length near 4.75 a 0 for most levels of treatment considered, and the r e values 
given in Table I refer to this local minimum. Explicit calculations show that these 
local minima are genuine and do not result from basis set superposition errors. 

The CC methods in Table I are listed in order of sophistication of the itera- 
tive component of the method. The CCSD method [18] includes all effects of the 
T\ and T 2 cluster operators, which include disconnected triple and quadruple ex- 
citation effects. The CCSD+T(CCSD) [7], CCSD(T) [ 8 ], CCSD+T*(CCSD) [2], 
CCSD+TQ(CCSD) [2], and CCSD+TQ*(CCSD) [2] methods extend CCSD by 
providing noniterative estimates of the effects of the triple excitation cluster oper- 
ator T 3 (ie connected triple excitation effects) and (for CCSD+TQ(CCSD) and 
CCSD+TQ*(CCSD) only) the effects of the quadruple excitation operator T 4 
(ie connected quadruple excitation effects) from CCSD T\ and T 2 amplitudes. 
All of these methods are correct through fourth order of perturbation theory. 
The CCSD(T) and CCSD+T*(CCSD) methods contain, in addition, nonitera- 
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tive estimates of different fifth-order terms, while the CCSD+TQ(CCSD) and 
CCSD+TQ*(CCSD) methods are correct through fifth order. The noniterative 
fourth- and fifth-order corrections to the CCSD energies are shown in Table II. For 
brevity, in the tables and the rest of the text we refer to the CCSD+T(CCSD), 
CCSD+T*(CCSD), CCSD+TQ(CCSD), and CCSD+TQ*(CCSD) methods as 
CCSD+T, CCSD+T*, CCSD+TQ and CCSD+TQ*, it being understood that 
the noniterative corrections axe evaluated with CCSD amplitudes. The various 
CCSDT-n methods [6,7,9,11] are all iterative schemes for estimating the effects of 
connected triple excitations; they differ by the number and type of terms included in 
the T z equation. The higher-numbered methods include more terms. The CCSDT 
method [9,10] includes all effects of the T\ % T 2 , and T z cluster operators, and is exact 
for a three-electron system. The CCSDT+Q(CCSDT) and CCSDT+Q*(CCSDT) 
methods [12] estimate effects of connected quadruple excitations (ie the T\ operator) 
noniteratively from converged CCSDT amplitudes. These estimates are also shown 
in Table II. These methods are hereafter referred to as CCSDT+Q and CCSDT+Q*, 
it being understood that the T 4 corrections are computed from CCSDT amplitudes. 

One obvious conclusion to be drawn from Table I is that, as we have observed 
in many other situations (see Ref. 19 and references therein), the multireference 
calculations reproduce the FCI results almost exactly. (We have not employed 
a multireference Davidson correction in the present calculations: as discussed in 
Ref. 19 such a correction generally overshoots the FCI result when only six electrons 
are correlated and a complete reference space is used, as in the present MRCI 
calculations.) The MRCI results themselves differ from FCI by tens of fiE h at all 
geometries, a discrepancy more than an order of magnitude smaller than for any 
of the other methods, and the geometry dependence of the difference is also much 
smaller. 

The CC methods that take some account of connected triple excitations offer 
substantial improvement over the CCSD method for Be 3 , all giving energies in much 
better agreement with FCI and all giving a local minimum which is absent from the 
CCSD potential. The noniterative CCSD+T, CCSD(T), and CCSD+T* methods 
are in very close agreement, because of the small values of the fifth-order terms 
E l g ] T and E l £ ] D (see Table II) which distinguish CCSD(T) and CCSD+T* from 
CCSD+T. These methods are also in very close agreement with CCSDT-1 and 
CCSDT-2, suggesting that iteration of the T 3 equation is not of great importance 
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in this case. CCSD(T) is slightly closer to FCI than is CCSD+T* since E ^ is 
much smaller than Ej)^ for this example (Table II). In other systems the reverse 
is sometimes true and the relative sizes of the two terms also vary with molecular 
geometry [2]. 

The more sophisticated noniterative methods CCSD+TQ and CCSD+TQ*, 
which are correct through fifth order, perform very well. They are in significantly 
better agreement with FCI than are the simpler noniterative methods and are com- 
petitive with CCSDT+Q and CCSDT+Q*, which are considerably more expensive. 
CCSD+TQ* is the more theoretically sound method for adding connected T 4 cor- 
rections to CCSD [2], and this method gives slightly better agreement with FCI 
than does CCSD+TQ. 

Although the CCSDT-2 and CCSDT-3 methods are theoretically more com- 
plete than CCSDT-1, they are in slightly poorer agreement with FCI. The reason 
for this is the presence of additional terms in the T3 equation which normally make 
positive contributions to the energy. This has been observed previously [9,11], and 
should not necessarily be viewed as a deficiency of these methods since the extra 
terms confer extra stability in some difficult cases. The CCSDT-4 method is in 
extremely good agreement with FCI, but this is fortuitous since the more com- 
plete CCSDT method is not in such good agreement. The CCSDT-4 method omits 
from the T 3 equation all nonlinear terms involving T3, and consequently tends to 
give energies slightly lower than CCSDT, and sometimes in better agreement with 
FCI [9]. 

The quadruples corrected CCSDT results are in slightly better agreement with 
FCI than are the CCSDT results, with the CCSDT+Q* method, like CCSD+TQ*, 
having an error of only a fraction of a m Eh. The T4 corrections are seen to be neg- 
ative, which is usually the case near equilibrium geometries. As discussed in some 
detail in Ref. 2, for correcting CCSDT energies the correction Q is to be preferred 
on mathematical grounds to Q*, the reverse of the situation when correcting CCSD. 
However, in the present case the Q* correction to CCSDT agrees better with FCI 
than does Q. 

The CC bond lengths given in Table I differ from the FCI value by more 
than 0.1 a 0 in some cases, which may appear disappointing in view of the perfect 
agreement between the FCI and MRCI values. This apparent deficiency should not 
be overemphasized, however, as the very flat potential around r e magnifies the effect 
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of small energy differences on the computed bond length. It is more reasonable to 
view all of the (triples-corrected) CC bond lengths as being in fair agreement with 
one another and with FCI — more importantly, the connected triples appear to 
be entirely responsible for the local minima in the CC potentials, as the CCSD 
potential has no local minimum at all. 

Additional insight into the importance of connected and disconnected higher 
excitations can be obtained from the Cl results given in Table I. The agreement 
between CCSDT+Q and CISDTQ is very close, and both agree well with FCI. 
However, disconnected quadruples obviously play a central role in achieving this 
agreement, as CISDT differs from FCI by more than any result in the table (other 
than CISD). From the difference between CISD and CCSD it appears that discon- 
nected triples and quadruples contribute some 13 m E\ to the energy: even in a case 
like Be 3 , in which near-degeneracy effects increase the importance of disconnected 
triples, it is likely that much of this difference is due to disconnected quadruples [20]. 
It is, of course, hardly feasible to include quadruple excitations in Cl calculations 
on most systems of interest because of the computational cost. 

In Table III we show the results obtained with the more realistic [4s 2p Id] basis 
set. The noniterative corrections to CCSD and CCSDT in this basis are shown in 
Table IV. With this larger basis set we do not have FCI values for comparison, but 
we may expect the MRCI method to produce closely similar results, and accordingly 
this is our reference point. The same general trends seen in Table I are observed 
here also. With the exception of CCSD and MBPT(4), the agreement between 
all methods is, on the whole, very good. The failure of MBPT(4) is a reflection 
of the need to include single and double excitation effects iteratively, which is not 
surprising for a system like Be 3 . The good agreement between the noniterative 
methods and their iterative counterparts shows, however, that it is not necessary 
to include connected triple and quadruple excitation effects iteratively. The full 
CCSDT, CCSDT+Q, CCSDT+Q*, CCSD+TQ, and CCSD+TQ* results are in 
excellent agreement with MRCI, although the (small) difference between MRCI 
and these CC methods has a relatively large variation with distance. 

In general, the variation among the computed bond lengths with different cor- 
relation treatments is smaller than is the case for the [4s 2p] basis set. Only the 
MBPT(4) value shows an error as large as 0.1 <io> an d most methods agree well with 
the MRCI result. 
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The binding energy results may be divided into two classes: apart from 

CCSD+TQ and CCSD+TQ*, the CC methods that account approximately for 
connected triples produce results of 12.0 kcal/mol or less, while CCSDT and 
CC methods that account for quadruple excitations yield binding energies around 
13.5 kcal/mol. The discrepancy between the various CC values and the MRCI result 
of 14.5 kcal/mol may seem laxge as a relative error, but comparisons in larger basis 
sets [4,5] show that the absolute error, at least at the CCSD(T) level, remains much 
the same, and hence in laxge basis sets the relative error in the binding energy is 
rather small. We stress that the [4s 2 p Id] binding energies are far from the basis set 
limit values: the MRCI result in a complete basis is estimated to be 24 kcal/mol [3]. 

We now analyze the apparent success of the simple noniterative CC methods 
CCSD+T and CCSD(T) in more detail. The performance of these methods is good. 
For the larger basis set, they remove about 87% of the error of the CCSD method, 
while the CCSDT method removes an additional 10%. Hence, triple excitation ef- 
fects beyond those included in CCSD(T) further increase the binding energy slightly, 
as do connected quadruple excitations, though not by as much. The fact that the 
residual triple and quadruple excitation effects are relatively small and of the same 
sign suggests that the good agreement of CCSD+T and CCSD(T) with MRCI and 
FCI is not fortuitous. It should be remembered, however, that the smallness of 
the residual triple excitation correction can be considered to be a consequence of 
a partial cancellation of several contributions to the energy of opposite sign (see 
Tables II and IV). To illustrate, we note that for the larger basis set the CCSDT 
method gives energies which are 1-1.7 m Eh. lower than the CCSD+T or CCSD(T) 
energies. However, using the CCSD+TQ and CCSD+TQ* estimates, the contri- 
bution to the energy of the E^ T term [2], arising from the presence of WT Z (see 
Ref. 2) in the T 3 equation, is about -3 m^, which is partially offset by the non- 
linear terms in the T 3 equation. From this analysis, therefore, one could argue that 
the closeness of the CCSD+T and CCSD(T) results to MRCI and FCI depends to 
some extent upon such cancellations. 

On balance, we believe that it is reasonable to apply the CCSD(T) method to 
larger Be clusters, as we have previously suggested [5]. Of course, in any application 
of the method to systems with substantial nondynamical correlation effects the 
potential for large residual triple and quadruple excitation contributions should 
always be borne in mind. We note that for Be 2 it is precisely the residual triple 
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excitation effects absent from CCSD(T) which are necessary in order to get a correct 
potential curve [1], although, as we have stressed above and elsewhere [3], there are 
few parallels between the bonding in Be2 and larger Be clusters. For Be3 the residual 
triple and quadruple excitation effects are small and of the same sign, but this will 
by no means always be the case. Other systems show quite different correlation 
effects [2,21]. 

Finally, considering that the CCSD+T method actually gives slightly better 
agreement with FCI and MRCI here, it seems tempting to use the CCSD-|- T method 
instead of CCSD(T). We do not recommend this since adding the positive fifth-order 
JjH. term will in general confer more stability, particularly in highly correlated sys- 

^ * . f5l 

terns, as has been discussed previously [2,4,8]. The cost of evaluating the E ST term 
is negligible compared with solving the CCSD equations or evaluating the fourth- 
order triples correction. CCSD+T* will have similar advantages over CCSD+T. 

IV. Conclusions 

Benchmark FCI, CASSCF/MRCI, and a variety of single-reference Cl and CC 
calculations have been performed on Be3 to determine the effects of connected 
triple and quadruple excitations and to compare with the simple noniterative CC 
method CCSD(T). The MRCI results are in excellent agreement with FCI. The most 
sophisticated single-reference Cl and CC methods, CISDTQ and CCSDT+Q, give 
very simil ar results and are in very good agreement with MRCI and FCI. It is found 
that connected triple excitations account for almost all of the difference between the 
CCSD and FCI energies, while the contribution of connected quadruple excitations 
is rather small, but also increases the binding energy. In this case, it appears not 
to be essential to include either T3 or T4 iteratively, since the noniterative methods 
give results very similar to their iterative counterparts. The failure of MBPT(4) 
shows, however, that the iterative inclusion of T\ and T 2 is important, at least as 
measured in the modest basis sets used in this work. The closeness of the CCSD(T) 
results to those of more extended treatments of electron correlation is seen not to 
arise from a cancellation of errors between neglected triple excitation contributions 
and higher excitations, in the sense that the residual triple and quadruple excitation 
corrections are both small and of the same sign. Our previous work [4] has also 
shown good agreement between CCSD(T) and extensive MRCI treatments for Be4, 
and it seems this also is a genuine agreement and not a cancellation of errors. 
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This reinforces our earlier optimism about the potential for success in applying the 
CCSD(T) method to larger beryllium clusters, although the possibility of terms 
neglected in this simple approach being significant should always be borne in mind 
when studying these highly correlated systems. 
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Table I. Comparison of SCF, MBPT(4), CC, Cl, MRCI, and FCI results for Be with the [4s 2 p] ANO basis set . 
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a Total energy in E*, differences relative to FCI (in parentheses) in mEfc, r and r e in a 0 . 
b Not computed. 

c No local minimum in potential. 


Table II. Fourth- and. fifth-order noniterative corrections to the CCSD and CCSDT 
energies at different bond lengths for the [4s 2 p] basis set 0 ’ 6 . 


Correction r =4.5 


4' 

-0.008055 

4S- 

-0.000071 

zpl 5 l 

&DT 

0.000472 

e [ t ] t 

-0.002334 

Q(T(CCSD)) 

-0.000403 

Q*(T(CCSD)) 

-0.000830 

Q(CCSDT) 

-0.000511 

Q*(CCSDT) 

-0.001133 


r =4.8 

r =5.1 

-0.006277 

-0.004579 

0.000015 

0.000031 

0.000261 

0.000109 

-0.001912 

-0.001495 

-0.000336 

-0.000253 

-0.000678 

-0.000501 

-0.000401 

-0.000305 

-0.000851 

-0.000625 


° Energies in Eh, r in a 0 . 

6 For definitions of the corrections see references 2 and 12. 
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Table III. Comparison of SCF, MBPT(4), CC, and MRCI results for Be with the [4s 2 p l<f] ANO basis set a . 
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Table IV. Fourth- and fifth-order noniterative corrections to the CCSD and CCSDT 
energies at different bond lengths for the [4s 2 p l<f] basis set®’ 6 . 


Correction 

r =4.2 

4 

-0.012680 

i?[ 5 i 

2j S t 

0.000090 

&DT 

0.000800 


-0.002973 

Q(T(CCSD)) 

-0.000251 

Q*(T(CCSD)) 

-0.000538 

Q(CCSDT) 

-0.000400 

Q*(CCSDT) 

-0.000905 


r =4.5 

r =4.8 

-0.011839 

-0.010345 

0.000167 

0.000234 

0.000658 

0.000493 

-0.002890 

-0.002706 

-0.000294 

-0.000163 

-0.000480 

-0.000351 

-0.000304 

-0.000206 

-0.000660 

-0.000444 


° Energies in Eh, r in a<>. 

i For definitions of the corrections see references 2 and 12. 
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